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We have calculated the plasmon modes in graphene double layer structures at finite temperatures, 
taking into account the inhomogeneity of the dielectric background of the system. The effective 
dielectric function is obtained from the solution of the Poisson equation of three-layer dielectric 
medium with the graphene sheets located at the interfaces, separating the different materials. Due to 
the momentum dispersion of the effective dielectric function, the intra- and inter-layer bare Coulomb 
interactions in the graphene double layer system acquires an additional momentum dependence-an 
effect that is of the order of the inter-layer interaction itself. We show that the energies of the 
in-phase and out-of-phase plasmon modes are determined largely by different values of the spatially 
dependent effective dielectric function. The effect of the dielectric inhomogeneity increases with 
temperature and even at high temperatures the energy shift induced by the dielectric inhomogeneity 
and temperature itself remains larger than the broadening of the plasmon energy dispersions due 
to the Landau damping. The obtained new features of the plasmon dispersions can be observed in 
frictional drag measurements and in inelastic light scattering and electron energy-loss spectroscopies. 

PACS numbers: 73.22.Pr; 73.20.Mf ; 73.21.Ac 



I. INTRODUCTION 

Graphene is a monolayer of carbon atomiP^with great 
potential for a new generation electronics^'. Charge car- 
riers in graphene are Dirac-like massless, chiral fermions 
that provide a unique two dimensional system with new 
many-body phenomengP that can critically influence the 
electronic properties'^ in graphene. An excellent tool 
for studying many-body interaction in graphene struc- 
tures are graphene double-layer systems (GDLS), re- 
cently realized in several experiments^^. 

Bare electron-electron interaction in graphene is de- 
scribed by the dimensionless fine-structure constant?^ 

a g = 7 — z ~ — (1) 
hvpe e 

with its value depending on the dielectric properties 
of the graphene surrounding environment via the ef- 
fective dielectric constant I (vf the electron velocity 
in graphene). In GDLS the two spatially separated 
graphene sheets are immersed in a nonhomogeneous three 
layered medium with background dielectric constants ei , 
62 and £3 of the contacting media, as shown in Fig. 1. 
In general, these dielectric permittivities differ substan- 
tially from each others in experiment. In an individual 
graphene sheet on top of a substrate with relative dielec- 
tric permittivity e^, an electron charge e behaves effec- 
tively as a charge with a value of 2e/(l + e^). This well 
known resultP^O has been applied in recent treatment of 
Coulomb drag in GDLS^, assuming that the effective 
permittivity in each graphene layer of the GDLS is given 
by the arithmetic average of its surrounding media, i.e. 
e 12 = (ei + e 2 )/2 and e 23 = (e 2 + £3)/2. In what fol- 
lows, we argue that such an approach in general is not 
applicable for the plasmon problem in GDLS. It neglects 
the momentum dispersion of the effective dielectric per- 



mittivity, e — e(q), an effect of the order of the inter- 
layer interaction. Its direct application would result in 
the disappearance of the linear dispersion of out-of-phase 
plasmon modes because the bare intra- and inter-layer 
Coulomb interactions remain not equal in the long wave- 
length limit. Furthermore, when an average spatially in- 
dependent dielectric permittivity is used^iHMl t describe 
GDLS, we find that it also provides an inadequate de- 
scription of the plasmon energy dispersions in GDLS. In 
the present paper we develop a consistent description of 
the plasmon modes in GDLS by making use of the ex- 
act solution of the Poisson equation for the electrostatic 
problem in three layer dielectric medium and by taking 
into account the momentum dispersion of the effective 
background dielectric function e(q) due to the finite thick- 
ness of the inter-layer barrier. We show that even in the 
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FIG. 1. A graphene double-layer system immersed in a three 
layered dielectric medium. The solid lines with spacing d rep- 
resent the graphene sheets 1 and 2 at the interfaces, separating 
different materials with the background dielectric permittivi- 
ties ei, £2, and £3. 

long wavelength limit, when the bare intra- and inter- 
layer Coulomb interactions remain equal, the dispersions 
of the in-phase and out-of-phase plasmon modes are de- 
termined by the largely different values of the spatially 
dependent effective dielectric permittivity. At finite mo- 
menta also the momentum dependence of the effective 
dielectric function e(q) results in significant changes of 
the dispersion relation of both plasmon modes which can 
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be probed in experiment. 



II. THEORETICAL CONCEPT 

We obtain the collective modes of plasmon 
excitations^ from the poles of the exact Coulomb 
Green function, V(q,uj). In double layer structures the 
bare Coulomb interaction, v(q), is a tensor with respect 
to the layer indices and in general represents three 
different interactions, the intra-layer, v±i(q) and V22(q), 
and the inter-layer, ^12(9) = ^21(9)- The kernel v(q) 
determines a standard matrix Dyson equation for the 
exact Coulomb Green function 



V(q, u) = v(q) + %) ■ n(<7, u) ■ V(q, u) 



(2) 



where Ii(q, uS) is the irreducible polarization function of 
the double layer electron system. It is seen from ^ 
that the poles of the exact Coulomb Green functions 
are given by the zeros of the scalar screening function 
e(q,w) = det|l — v(q) ■ II(g,a;)|. In GDLS even for the 
inter-layer spacing of 2-3 nm, the Dirac carriers in the 
graphene sheets are coupled via inter-layer Coulomb in- 
teraction and the tunneling between the layers is insignif- 



icant. Hence we can neglect the non-diagonal elements 
of the polarizability tl(q,u>). It is also sufficient to re- 
strict ourselves to consider the polarizability within the 
random phase approximation where flu (q, uj) — II^ (q, lu) 
and II22 (<jf, w) = fl^g, a>) and the noninteracting polar- 
ization functions II° 2 (<7,u;) are given by the bubble dia- 
grams in the respective graphene sheets. This approach is 
well justified for weakly interacting GDLS and its screen- 
ing function can be represented as 

e(q,u>) = e 1 (q,uj)e 2 (q,uj) - wi 2 (g) 2 n°((j, uj)U° 2 (q, u>) (3) 

where £1,2(9, w) = 1 — Uii,22(<?)n? 2(91^) are the screen- 
ing functions in each graphene layer. Despite the external 
similarities, the screening function of GDLS differs essen- 
tially from that of the usual two-dimensional electron gas 
in semiconductor nanostructures. In addition to the new 
properties of H°(q,ui) arising from the unique Dirac-like 
energy bandstructure and the chiral nature of the mass- 
less carriers, we show here that the heterogeneity of the 
dielectric background in GDLS plays an important role in 
determining the many-body Coulomb interaction effects. 

Direct calculations of the electrostatic problem in the 
dielectric environment consisting of three contacting me- 
dia with different dielectric constants e 1; e 2 , and e 3 give 
the following formula for the Coulomb potential 



ip(z,z ) 



2ne 2 



-q\z-z \ 



q^2 

(^2 - £3) (ei sinhqzo + e 2 coshqz ) e ? ( z ~ d ) — (ei — e 2 ) (e 2 coshq(d - z ) + £3 sinhq(rf — z )) e~ 



£2 (ei 



£3) coshdg + (£163 -I- e|) sinhgd 



(4) 



where z and z are, respectively, the positions of the 
source and test charges. The bare intra-layer interactions 
in graphene layers 1 and 2, located at the two interfaces 
with z = and z = d, are given by vn(q) = eip(0, 0) 
and t>22 = e-fid, d) wh ile the bare inter-layer interaction 
V12 (q) = ec^(0, dp^l 2 !!. Hence, we arrive at the following 
three effective dielectric functions in the GDLS 



(5) 
(6) 
(7) 
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that determine the strength of the intra- and inter-layer 
bare Coulomb interactions in momentum space 



Vij(q,d) 



2Tre 2 
ql VJ (qd) 



(8) 



Here i,j = 1, 2 are the graphene layer indices. It is seen 
that in the long wavelength limit all three interactions 



are determined by the same effective dielectric constant, 
given by the arithmetic average of the top and bottom 
surrounding media in GDLS, £13 = (ei + £3)/2, and does 
not depend on the dielectric constant, €2, of the middle 
medium. 

Further we exploit the above formulae for bare inter- 
action Vij(q) to find the plasmon spectrum in GDLS. As- 
suming for simplicity that the density is balanced in the 
GDLS, we can rewrite the screening function ^ as 

e(q,Lj) =a(ll (q,u;) -IL+{q,d)) (lL°(q, u ) - n"( g ))(9) 

where a(q, d) — Vn(q, d)v2 2 {q, d) — v\ 2 {q, d) and the aux- 
iliary functions 



Here for brevity we omit the arguments of the bare in- 
teractions. In the limit of vanishing q we have 



V+{q,d) = ^q, lT{q,d) 



Aire 2 



2ne 2 d 



(11) 
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FIG. 2. The effect of nonhomogenous dielectric background on the plasmon dispersions in GDLS at low temperatures T = O.lTp. 
The upper (lower) curves represent the optical (acoustical) plasmon branches with the square-root (linear) dispersions. The 
solid curves are calculated for the plasmon modes in GDLS with nonhomogenous dielectric background, consisting of three 
layers of Si02 (esi0 2 — 3.8), AI2O3 (eAi 2 o 3 = 6), and air. The dashed curves correspond to the plasmons in GDLS with 
homogenous dielectric background with an average permittivity e = 2.4. In the left, mid, and right figures we plot, respectively, 
the plasmon energy dispersions for three values of the inter-layer spacing d = 5 nm, 10 nm, and 30 nm. The doping level 
corresponds to the carrier densities m = 71,2 = 10 12 cm -2 . The thick dot-dashed lines show the boundaries of inter- and 
intra-chirality particle-hole continua where plasmons are Landau damped. The shaded areas in the mid plot represent the 
broadening of the energy dispersions of respective plasmon modes. 



In unbalanced GDLS instead of the above formulas, one 
can find 



n+(a) = 



£13 



<1 



2ne 2 yfn\ 



n-(d) = 



C-2 



2ire 2 2dy/niU2 

and the dielectric inhomogeneity has a similar effect also 
in unbalanced systems. Here the partial dimensionless 
densities are defined as n ii2 = Ni^/N with N lt2 the 
carrier densities in each graphcnc layers and N = N1+N2 
stands for the total density in GDLS. 

Making use of the expression (8) from Ref. [55] for the 
zero temperature exact Lindhard polarization function in 
graphene, we find from the zeros of the screening function 
in Eq. ([9| the optical and acoustical plasmon modes in 
GDLS with the following square-root and linear energy 
dispersions 



/ ge 2 v F k F 



£13 



9 , u-(q) 



1 + q FF d 

V 7 ! + 2q TF d 



v F q . (13) 



Here g accounts for the spin and valley degeneracy in 
graphene and qr F = ge 2 k F /(e2E F ) is the Thomas-Fermi 
screening wave vector in graphene with k F defined by the 
single layer density. These plasmon dispersions, derived 
for GDLS in Ref. f2T)l f see the Note added), have the same 
form as for semiconductor two dimensional systems^. It 
is seen, however, that in contrast to the previous treat- 
ments, adopting a model of homogenous background di- 
electric environment for the GDL£pUH24j, ^ e pi asmon en _ 
ergies ( 13 1 of the optical and acoustical modes in GDLS 



are determined by the different background dielectric 
permittivities, which is a direct consequence of the be- 



havior of the auxiliary functions from (11) in the long 
wavelength limit. 



III. RESULTS AND DISCUSSIONS 

It follows from the above formulas that the energy of 
the in-phase optical plasmons is given by the arithmetic 
average £13 of the dielectric constants of top and bot- 
tom dielectric media in the GDLS and is independent 
of the middle layer permittivity, e 2 , while the energy of 
the out-of-phase acoustical plasmon modes is indepen- 
dent of ei, £3 and depends only on the dielectric con- 
stant e 2 . In recent experimental samples of GDLS 15 " 17 
the top graphene sheet is surrounded by air while the 
inter-layer barrier has a relatively large dielectric per- 
mittivity. In samples where the bottom graphene layer 
in GDLS lies on a boron nitride substrate, used in the 
experimental setup in Ref. [TTl ei = e 2 = £bn ~ 5 and we 
have £13 w 3, which is by a factor 1.7 smaller than e 2 = 5. 
In samples on a Si02 substrate with moderate dielectric 
effects, studied in Ref. [Pol and [TBI ^ = esio 2 ~ 3.8 and 
(2 = eAi 2 o 3 ~ 6 and we have £13 ss 2.4, which is smaller 
than £2 = 6 by a factor of about 2.5. 

Thus, the inhomogeneity of the background dielectric 
environment in GDLS can significantly alter the disper- 
sions of plasmon modes in the long wavelength limit by 
reducing essentially the energy of the out-of-phase plas- 
mon branch (see numerical calculations below). In par- 
ticular, this can lead to strong modifications of the tem- 
perature dep enden ce of the plasmon-mediated Coulomb 
drag in GDLS^211tlIl At low temperature the plasmon- 
mediated drag is mainly determined by the acousti- 
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FIG. 3. (left) Temperature effect on the in-phase and out-of-phase plasmon modes in GDLS with nonhomogenous dielectric 
background. The dashed and solid lines correspond to the T = O.ITf and T = T F . (mid) and (right) The effect of nonho- 
mogenous dielectric background on the plasmon dispersions in GDLS at high temperature T — T F . The plasmon dispersions 
are shown (mid) in the exactly balanced, ni = ni = 10 12 cm -2 , and (right) in the completely unbalanced, «2 = 0, GDLS. In 
all figures d — 10 nm and other parameters and notations correspond to Fig. [2] The shaded areas in the mid plot represent the 
broadening of the energy dispersions of respective plasmon modes. 



cal plasmon modes, which will be easy to excite ther- 
mally because they become closer to the boundary of the 
particle-hole continuum. As a result the upturn tempera- 
ture of the plasmon enhancement will be shifted to lower 
temperatures with the higher enhancement peak of the 
drag resistivity. 

The momentum dispersion that the effective dielectric 
functions Cij(q) exhibit in Eqs. (5j7) modifies substan- 
tially also the double-layer plasmon dispersions at finite 
values of q. In Figs. [2] and [3] we calculate the plasmon 
dispersions at two different temperatures T = O.ITf and 
T = Tp in the range of momenta, < q < 2kp, from 
the zeros of the screening function ^ by making use 
of the exact semi-analytical formulas from Ref. [5H for 
the finite temperature polarization function H^(q,u\T) 
of graphene. 

As seen in Fig. [2] the account for the inhomogeneity of 
the dielectric background in GDLS suppresses strongly 
the out-of-phase acoustical plasmon mode. The largest 
energy difference in comparison with the homogenous 
case is achieved nearly at the boundary of the inter- 
chirality particle-hole continuum for q ~ k F . The effect 
of inhomogeneity of the dielectric background on the en- 
ergy of in-phase plasmon modes starts to be significant 
at finite values of q ~ 0.2kp. It increases with q so that 
the dispersion curve remains approximately parallel to 
the boundary of the intra-chirality subband particle-hole 
continuum at momenta q ~ 2kp. In contrast, the energy 
deviation of the out-of-phase plasmon modes due to the 
dielectric inhomogeneity shows nonmonotonic behavior 
and for large values of q ~ kp starts to decrease with q. 
It is seen also that the dispersion curves of the in-phase 
and out-of-phase plasmon modes exhibit an usual behav- 
ior and with an increase of the inter-layer spacing two 
modes become degenerate outside the particle-hole con- 
tinua. The shaded areas in the middle figures of Figs. [2] 



and [3] represent the broadening of the respective plasmon 
dispersions in the complex uj plane. Instead of the abso- 
lute value of 3II (g,w), the broadening T±(q) is the real 
measure of the plasmon Landau damping and is given 
bjESl 



r±(g) 



314° (g,w) 
dm°(q,uj)/duj 



(14) 



Far from the boundaries of the particle-hole continua, the 
plasmon modes are strongly Landau damped and the dis- 
persion curves do not represent well defined elementary 
excitations. With an increase of T the particle-hole con- 
tinua "degrade" while the damping of plasmon modes 
with smaller values of q increases. It is seen, however, 
that even at high temperatures T — Tp the broadening 
of the plasmon dispersions for q < kp is smaller than the 
shift of the energy dispersions, induced by the combined 
action of the dielectric background inhomogeneity and 
temperature. This is true both for the optical and acous- 
tical plasmon modes. In this small q and small u) region 
QTI (a;, q) is small and because its contribution to the 
real part of the screening function in Eq. (|3| is quadratic, 
the plasmon dispersion relations obtained from the zeros 
of the factors of the screening function in (|9| give the 
same result against the background of the broadening of 
plasmon dispersions, linearly varying with £sll (<7, id). 

In Figs. [3] we study the combined effect of finite tem- 
peratures and the inhomogeneous dielectric background 
on the energy dispersions of double-layer plasmon modes. 
It is seen in Fig. [3jleft) that the plasmon energy of both 
modes increases significantly with T. The comparison 
of Figs. J2^mid) and [3|mid) shows that the effect of the 
dielectric background inhomogeneity becomes stronger 
with T. As seen from Figs. 3][mid) and fright) at fi- 
nite temperatures T — Tp the plasmon modes behave 
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almost in the same way in two extreme regimes of the 
completely balanced, n± = ri2, and the completely un- 
balanced, n 2 = 0, GDLS. This is in stark contrast to 
the behavior that plasmons exhibit at T = At the 
same time Fig. fright) shows that the effect of dielectric 
background inhomogeneity in the unbalanced GDLS is 
similarly strong as in the balanced GDLS. 



IV. SUMMARY 

We study the effect of the dielectric background in- 
homogeneity on the plasmon modes in GDLS at finite 
temperatures. It is found that using a spatially aver- 
aged dielectric permittivity to describe dielectric proper- 
ties of GDLS provides an inadequate description of the 
plasmon energy dispersions. We obtain the effective di- 
electric functions of GDLS from the exact solution of the 
Poisson equation for a three-layer dielectric medium with 
the graphene sheets located at the interfaces, separat- 
ing different materials. It is shown that the momentum 
dispersion of the effective dielectric function results in 
an additional momentum dependence of the intra- and 
inter-layer bare Coulomb interactions in GDLS-an effect 
of the order of inter-layer interaction itself. As a result 
the obtained dispersions of the in-phase and out-of-phase 
plasmon modes are largely determined by the different 
values of the spatially dependent effective dielectric per- 
mittivity. 

Here we use the temperature dependent exact Lind- 
hard polarization function, which allows us to find the 
plasmon dispersions at finite temperatures. Our calcula- 
tions show that the effect of dielectric background inho- 
mogeneity increases with temperature and even at high 
temperatures, T ~ Tp, the energy shift induced by the 
dielectric inhomogeneity and temperature itself remains 
larger than the broadening of the plasmon energy dis- 
persions due to the Landau damping at finite tempera- 
tures. This combined effect of finite temperatures and 
the dielectric background inhomogeneity has been dis- 
cussed thoroughly in our paper by comparing carefully 
the dispersion relations, calculated for graphene dou- 
ble layers embedded in a three-layer nonhomogeneous 
dielectric medium with that obtained for the homoge- 
neous one. We find a strong effect of temperatures on 
the acoustical plasmon mode in a completely unbalanced 
system-the acoustical mode becomes separated from the 
top of the electron-hole continuum and this favors pos- 
sibility to observe it. Our numerical calculations have 
been carried out for realistic experimental samples on a 



Si substrate with three different dielectric constants for 
the background nonhomogeneous medium for interlayer 
spacing d — 5, 10, and 30 nm. The effect of dielectric in- 
homogeneity is found to be the largest for samples with 
d = 10 nm at finite momenta q < kp. 

The predicted new features of the double-layer 
plasmon dispersions can be observed in high tem- 
perature measurements of frictional drag in GDLS 
and in such plasmon experiments as plasmon- 
enhanced p hotolu minescence 3 -^, electron energ y-loss 
spectroscopies^, inelastic light scattering } 37 ! 38 ! mea- 
surements. 

Note added — During the preparation of our 
manuscript when its main results have been already 
summarized in the abstract submitted to the March 
MeetingpSl, we became aware of two other preprints 
(Refs. l29l and l40|) on the double-layer plasmon modes in 
graphene structures. In the present paper we give par- 
ticular emphasis to the effect of inhomogeneity of the di- 
electric background on the plasmon dispersions in realis- 
tic GDLS and use a different formalism from that used in 
Refs. [inland HOI which allows us to treat additionally the 
effect of finite temperatures on the energy dispersions of 
the plasmon modes. The main focus of Ref . [2"51 has been 
the subtle point related to the use of the exact Lindhard 
polarization function in obtaining the plasmon velocity 
in the long wavelength limit. The authors derived first 
the analytical formula for the acoustical plasmon veloc- 
ity thereby correcting the previously reported result from 
Ref. [5T] The effect of dielectric inhomogeneity, which is 
due to the renormalization of the bare Coulomb interac- 
tions, in Ref. [3nihas been discussed mainly in connection 
with the effect of the locked acoustical mode on top of 
the particle-hole continuum in topological insulators. We 
find, however, that the numerical calculations presented 
in Figs. 2-4 of Ref.l29lfor the interlayer spacing d = 0.335 
nm are not adequate. For such a small value of the barrier 
thickness the tunneling between graphene layers cannot 
be neglected and the Eq. (8) in Ref. [25] is inapplicable. 
Although the authors of Ref. 40! have reported the main 
formulas for the three layer dielectric medium, their main 
focus has been the study of the near field amplification 
in GDLS with a homogenous dielectric background. 
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